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We put forward a scheme to study the anisotropic magnetic couplings in SrilrCL by mapping fully relativistic 
constrained noncollinear density functional theory including an on-site Hubbard U correction onto a general 
spin model Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, 
spin and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy 
and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in 
Sr 2 Ir 04 arises from the interplay between structural distortions and the competition between isotropic exchange 
and DM interactions. A complete magnetic phase diagram with respect to the tetragonal distortion and the 
rotation of IrC>6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic 
transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal 
weak ferromagnetic moment below a critical octahedral rotation. 


In weak ferromagnetic materials the subtle interplay among 
different types of magnetic interactions can cause the forma¬ 
tion of complex canted spin structures involving the so-called 
Dzyaloshinskii-Moriya (DM) effect, arising from the coupling 
between the spin and orbital angular momenta Hi. Intense 
research was done in this field in the last few years, motivated 
by the foreseeable applications in storage technology and by 
the air of mystery enveloping the quantum-mechanical origin 
of DM structures M- A crucial aspect of the DM systems 
is the entanglement between structural distortions and mag¬ 
netism, which could be exploited as a way to tune the spin 
texture by modifying the structure upon external stimuli such 
as pressure and strain 001. 


The cross coupling between the different electronic, lattice 
and spin degrees of freedom is particulary rich in iridates. 
Here, the spin-orbit coupling (SOC), electron-electron corre¬ 
lations, and spin-exchange interactions operate with compara¬ 
ble strengths and gives rise to a large variety of exotic states 
SE1 . The most striking example of this class of materials is 
the layered perovskite SmlrCL, characterized by a novel rel¬ 
ativistic Mott insulating state EMI El! [l7ll and an unusual 
in-plane canted antiferromagnetism JA.FM) with a weak net 
ferromagnetic (FM) component |4, 181. The small electronic 
gap (=s 0.3 eV 117]) is opened by modest Hubbard interac¬ 
tions (U ~ 1.5-2 eV [IS]) and by the strong spin-orbit cou¬ 
pling (% soc =» 0.5 eV 12011 ) which effectively narrows the d 
orbital bandwidth and give rise to an ideal 7 e ff= 1/2-like state 
im imm. This is considered to be robust despite the pres¬ 
ence of noncubic structural distortions I 23 I 1 . Neutron diffrac¬ 
tion experiments indicate that the lrO f , octahedra are rotated 
by a =11.5° and elongated in the c direction {c/a ~ 1.04) 
[24], generating the enlarged \[2ax2c I4\/acd tetragonal cell 
shown in Fig. [O' a). The spins, coupled with the orbital mo¬ 
ment, follow the rotation of the octahedra and establish an in¬ 
sulating canted AFM ordering with a canting angle 6 ~ 12.2° 


[25] 1261] [Figs- Eta) andQJb)]. The formation of the canted 
AFM ordering is believed to arise from the entanglement of 
structural distortions and SOC, which affects the balance be¬ 
tween the DM anisotropic interaction and the isotropic mag¬ 
netic coupling. A simultaneous treatment of all - lattice, spin 
and orbital - effects is therefore essential for an accurate ac¬ 
count of magnetic interactions in this type of magnetically 
canted structures. 


The aim of our study is to propose a general scheme to treat 
and explain canted magnetic states in spin-orbit coupled sys¬ 
tems fully ab initio. We do this by introducing a series of 
directionally constrained non-collinear spin setups and self- 
consistent total energy calculations. We determine the magni¬ 
tude of all relevant magnetic interactions in S^IrC^, including 
the isotropic exchange, the single-ion anisotropy and the DM 
coupling, and study the influence of structural perturbations 
(tetragonal distortion and octahedral rotation) on the spin or¬ 
dering. 


The theoretical modelling of canted states can be conducted 
by model Hamiltonian or first-principles approaches. Based 
on a model Hamiltonian approach lackeli and Khaliullin de¬ 
termined that magnetic interactions in SiMrCL are governed 
by the lattice geometry enabling the possibility to tune the 
spin ordering by small structural perturbations. For instance, 
they found that a spin-flop transition to z-collinear AFM or¬ 
dering can be obtained above a critical tetragonal distortion 
strength 0. The calculation of the relative strength between 
the isotropic and anisotropic interactions is theoretically and 
computationally challenging. While the isotropic contribu¬ 
tions can be estimated accurately from first-principles within 
density functional theory (DFT), beyond-DFT methods and 
multiple scattering theory 127,128}], the accurate estimation of 
the DM coupling is extremely difficult. Using a phenomeno¬ 
logical microscopic model with a set of optimally chosen pa¬ 
rameters Kim and colleagues found a large DM coupling pa- 
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rameter D — 26.2 meV, fingerprint of a strong SOC, and a \D\ 
/ |/| ratio of 0.34, yielding 9 ~ 9.3° in good agreement with 
experiment 0. The ab inito calculation of D for periodic 
systems can be achieved by perturbation theory and Green’s 
function technique 110, scattering theory |3J |. or within 
constrained non-collinear spin DFT 1 51 33tl . through the spin- 
orbit induced corrections to the energy of spiral spin-density 

waves SI- 



FIG. 1. (color online) (a) The unit cell of S^IrCh containing 56 
atoms, (b) Projection onto the ab plane showing the octahedral tilting 
angle a and the canting angle 6 (c) Metal (M) to insulator (I) tran¬ 
sition phase diagram as a function of the Hubbard U and the SOC 
strength O. Filled circles are the computed (U, fl) points. The right 
(coloured) bar indicate the value of the gap (in eV) in the insulating 
regime. 

Here we follow an alternative route based on mapping 
self-consistent DFT total energies of several non-collinear 
magnetic configurations, obtained by constraining the mag¬ 
netic moment along specific directions, onto a general spin- 
dependent Hamiltonian. In our approach, we consider the 
cross-coupling between spin, structural relaxations and spin- 
orbit interaction fully self consistently and within the same 
theoretical framework. We clarify that the origin of the canted 
AFM state in SroIrC^ arises naturally from the fine balance be¬ 
tween the isotropic ( JS 2 - -0.32 meV), and DM ( D : S 2 = -0.25 
meV) couplings, and construct a detailed magnetic phase dia¬ 


gram showing the transition from the the canted ground state 
to different types of collinear magnetic states as a function of 
the tetragonal distortion c/a and the octahedral rotation angle 
a entirely ab initio. 

We use the Vienna Ab initio Simulation Package (VASP) 

135, 36] and adopt a fully relativistic (i.e., including SOC) 
DFT scheme within the gradient corrected approximation 
(GGA) with the inclusion of an on-site Hubbard U correc¬ 
tion [37] to account for electron screened exchange in the Ir 
d manifold. The reciprocal space is sampled with a 7x7x3 k- 
mesh and energy cut off up to 800 eV with a high convergence 
precision of 10 8 eV is necessary to achieved well converged 
results in the meV scale. The atomic positions are optimized 
with the lattice parameters fixed at the experimental values 
10 - 

This GGA+C/+SOC approach correctly describes the 
canted insulating ground state of Srilr 04 in line with previous 
LDA+f/+SOC studies fioh . The gap is opened by the inter¬ 
woven action of SOC and U as shown in the SOC /U metal- 
insulator phase diagram of FigHJc), constructed by tuning the 
SOC strength through the scaling factor Q (£2 = 1 correspond 
to the self-consistent SOC strength) for values of U ranging 
from 0 to 2.5 eV. From the metal-insulator phase diagram it 
can be seen that neglecting either the SOC (Q = 0) or U (U — 
0) Sr 2 lr 04 remains metallic: SOC and U alone cannot activate 
any metal-to-insulator transition (MIT). The MIT is driven by 
the cooperative action of SOC and U. As a general tendency 
we observe that the stronger the SOC strength, the smaller is 
the critical interaction U required for opening the gap. By 
fitting the U to the experimental bandgap (~ 0.3 eV H) 
for £1 = 1 we obtain an optimal U of 1.6 eV, similar to the 
constrained random phase approximation estimation, 1.96 eV 
0 - From the total energy difference with and without the 
inclusion of SOC we have estimated a SOC energy % soc of 0.7 
eV/Ir, in line with electron spin resonance measurements lf38j| . 
which is only about half of the computed value of U. These re¬ 
sults provide clear evidence for the relativistic Mott-Hubbard 
character of the insulating state in Sr 2 lr 04 . Moreover, the 
strong ligand field energy 10 Dq (3.7 eV, computed from the 
splitting between the e g and ? 2 g orbitals) pushes down the 
7 e ff=3/2 manifold and promotes a J,.n = 1/2 like state compat¬ 


ible with previous interpretations lllOll . Our estimated values 
of U, % soc and 10 Dq are collected in Table. Q] along with avail¬ 
able literature data. 

Moreover, our approach successfully reproduces the in¬ 
plane canted AFM state with a weak net moment of about 
0.08 1-IbI Ir, in good agreement with the experimental values 
of 0.06~0.1 jUfi/lr [l2, JJ, 0. The components of the or¬ 


bital (0.27 nsl Ir) and spin (0.12 /r^/Ir) moments are antiferro¬ 
magnetic ally aligned along the b axis and ferromagnetically 
canted with respect to the a axis [see Fig. |T]b)] . The cal¬ 
culated total canting angle is 8 - 12.3° (with similar spin, 
9 S = 14.4°, and orbital, 6 0 = 11.2°, components) is in ex¬ 
cellent agreement with the experimental value of 12.2 ° 0. 

After providing a detailed description of the ground state 
properties of Sr 2 lr 04 , we can now address and interpret the 
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TABLE I. Hubbard U (eV), ligand field energy 10 Dq (eV), spin- 
orbit coupling energy % soc (eV), isotropic AFM exchange JS 2 (meV), 
single-ion anisotropy K (meV), and the DM exchange D-S 2 (meV) 
of Sr 2 lr 04 . Available literature data are given with respect to the 
calculated S =0.12 p B . 
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origin of this exotic anisotropic magnetic structure by an ab 
initio evaluation of the magnetic coupling parameters within 
the framework of constrained DFT+f/ for noncollinear mag¬ 
netism. To this end, we consider a penalty contribution to the 
total energy which fixes the local moment into a specific di¬ 
rection 


E = £o + X W - M°(M° • M,)] 2 , (1) 

i 

where E a is the usual DFT energy without any constraint and 
the second term represents the penalty energy due to the non¬ 
collinear directional constraint. M° is a unit vector along the 
desired direction of the magnetic moment at site i and M, is 
the integrated magnetic moment inside the Wigner-Seitz cell 
around atom i, whereas the parameter A controls the penalty 
energy contribution, which becomes vanishingly small by in¬ 
creasing A. Thus, by progressively increasing A, one con¬ 
verges towards the DFT total energy for a given constrained 
noncollinear magnetic configuration. We found that T=10 is 
generally sufficient to reduce the penalty energy to < 10 5 eV. 

By varying the spin canting angle 9 S within the ab plane 
from 0 (afe-collinear AFM) to 45° (perpendicular alignment 
of nearest neighbor spin moments) we obtain the DFT total 
energy curve shown in Fig. [2] To calculate the magnetic cou¬ 
pling parameters, we start from the conventional expression 
for a spin-dependent energy 

au = - Z j '/ s ' ■ s j + Z 4, ( s.o + Z D 0 ■ t s ‘ x S A ( 2 ) 

i< j i i< j 

where the first, second and last term represent the 
isotropic AFM exchange, the single-ion anisotropy (SIA) and 
Dzyaloshinskii-Moriya interaction, respectively. Considering 
that Sr 2 lrC >4 adopts a quasi two-dimensional geometry (see 
Fig ID- the interlayer hopping is much smaller than the in¬ 
plane one and can be neglected; we therefore consider the 
in-plane nearest-neighbor hopping only. Also, within the in¬ 
plane canted AFM state the DM anisotropy has the z compo¬ 
nent only. In addition, the single-ion anisotropy reduces to 
e an (S) — Kcos(49 s ) as a result of the tetragonal symmetry 


in the lattice. After taking the sum over all the Ir ions in the 
conventional unit cell, Eq. (0) finally reduces to, 

A E = 16JS 2 cos{29 s ) + 8Kcos(49 s ) + 16D Z S 2 sin(20 J ) (3) 

where J, K and D z are the nearest-neighbor isotropic ex¬ 
change, SIA and DM coupling parameters, respectively, and 
S is the magnitude of the in-plane spin moment (see Fig. Q}. 

The magnetic coupling parameters J, K and D z can be ex¬ 
tracted by fitting the first-principles calculated data using the 
model Eq. (0). The obtained results, JS 2 = -0.32 meV, K 
— -0.1 meV, and D Z S 2 — -0.25 meV, collected in Table Q] 
compare well with phenomenological data available in liter¬ 
ature and elucidate the canted state completely ab initio. The 
AFM isotropic exchange, favoring a collinear alignment of the 
spins, is the dominating interaction and aids the stabilization 
of the AFM ordering. The DM interaction, which prefers an 
orthogonal coupling of the spins, is only slightly weaker than 
the isotropic exchange and assists the formation of a canted 
spin arrangement. The negative sign of D z indicates that the 
D vector is antiparallel to the vector product [S, x Sj] in order 
to reduce the energy of the system. The SIA is smaller than the 
isotropic and DM interactions, but still favors in-plane mag¬ 
netism at small canting angles, as can be seen from the de¬ 
composition of the energy curve into its isotropic, DM, and 
SIA contributions (Fig. 0. As expected, the evolution of the 
DM and isotropic energies as a function of 6 S follows a dif¬ 
ferent trend. The DM term become stronger by increasing 
9 S with a peak at 45° (perpendicular spin ordering) whereas 
the AFM isotropic exchange is more active for small canting 
angles and reaches its maximum strength at the onset of the 
collinear regime ( 6 S = 0°). The formation of the canted in¬ 
plane state at 9 S = 14.4° is the results of the subtle competition 
between these three terms, JS 2 , D Z S 2 and K. 



Spin canting angle 0 (deg.) 


FIG. 2. (color online) The spin-dependent total energy as a function 
of the in-plane spin canting angle 0, decomposed over the isotropic 
AFM exchange, SIA and DM contributions. 










4 


Having identified the origin of the anisotropic magnetism, 
we are now in a position to define the role of the structural dis¬ 
tortions on the magnetism and in particular on the way the c/a 
ratio and the IrC>6 octahedra angle a mediate the spin canting. 
In fact, due to the strong SOC, structural distortions inevitably 
affect the shape of the t 2 g orbitals and directly influence the 
magnetic interactions. 

Figure [3] shows the evolution of the canting angle 9 as a 
function of the tetragonal distortion (controlled by the octa¬ 
hedral c/a) and octahedral rotation angle a. The data, de¬ 
composed over the spin, orbital and total contributions, are 
obtained by varying c/a and keeping a fixed to its optimized 
value (11.5°) and vice versa. The plots show a strong coupling 
between the lattice and magnetic channels and very similar 
curves are found for the spin and orbital part. 

For a fixed a, the canting angle 9 exhibits a parabolic be¬ 
haviour as a function c/a [Fig [2a)]. At the cubic limit of 
c/a = 1, 9 (w 5.1°) does not vanish, indicating that the strength 
of the DM coupling is strong enough to produce magnetic 
canting in the presence of a finite octahedral rotation angle. 
For larger tetragonal distortions c/a > 1.09, a spin-flop non- 
collinear to c-collinear transition occurs, supporting the mi¬ 
croscopic predictions by Jackeli and Khaliullin Hi . 

When the c/a is fixed [Fig. [3fb)], the magnetic canting an¬ 
gle 9 follows the octahedral rotation angle a due to the strong 
SOC, confirming the robustness of the 1/2-like state even 
in the presence of a tetragonal distortion, in agreement with 
experimental observations [26|]. For a smaller than the criti¬ 
cal value of 4 °, the weak ferromagnetic moment is quenched 
and Sr 2 lr 04 converts to an in-plane collinear AFM insulator, 
consistent with the pressure-induced canted-to-collinear tran¬ 
sition inferred from the high-pressure experiments |@|. 

By conducting a series of additional calculations for several 
c/a and a values we obtained the ab initio magnetic phase dia¬ 
gram shown in Fig. [4] illustrating the regions of stability in the 
c/a-a space of the three different magnetic phases introduced 
in Fig. [3 ab-c anted AFM, a/?-collinear AFM and c-collinear 
AFM. Our results demonstrate that structural modifications 
can drive two types of magnetic transitions: (i) a spin-flop 
Morin-like transition, which is controlled by the tetragonal 
distortion and occurs at any value of a, and (ii) an intraplane 
canted to collinear transition emerging only at small roations 
(.a < 7° and around the cubic phase (c/a = 1 ± 0.06). The 
geometrical tunability of the magnetic ordering in Sr 2 lr 04 is 
made possible by the coexistence and competition of AFM- 
isotropic and DM interactions, the latter being particulary 
strong due to the large spin-orbit strength. 

In conclusion, by combining relativistic noncollinear con¬ 
strained DFT+f/ calculations with the spin-dependent energy 
model, we have evaluated the isotropic and Dzyaloshinskii- 
Moriya exchange parameters of SiMrC^ and built a struc¬ 
ture/magnetism phase diagram fully ab initio, involving two 
types of structure-induced canted to collinear magnetic tran¬ 
sitions. The magnetic interaction parameters are strongly af¬ 
fected by the lattice distortions, appealing for the direct ma¬ 
nipulation of complex magnetic states by small structural per¬ 


turbations. 
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FIG. 3. (color online) Change of the magnetic canting angle 9 (spin, orbital and total moments) as a function of the structural distortions, (a) 6 
vs. c/a at the optimized octahedra rotation angle a = 11.5°; (b) IrOj angle a vs. 6 at the optimized tetragonal distortion c/a = 1.04. Shaded 
areas indicate regions with a different spin ordering (nb-canted, ah-collinear and c-collinear, schematically represented in the insets). 



FIG. 4. (color online) Ab initio magnetic phase diagram of SrilrCL 
in the c/a-a space. Open (blue) circles refer to points selected from 
the calculated data to demarcate the phase transition boundary. The 
insets illustrate the specific magnetic ordering in the three regions. 
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